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Abstract 

A simple approximate solution for the quantum-mechanical quartic oscillator V = 
m 2 x 2 + gx A in the double-well regime m? < at arbitrary g > is presented. It 
is based on a combining of perturbation theory near true minima of the potential, 
semi-classical approximation at large distances and a description of tunneling under 
the barrier. It provides 9-10 significant digits in energies and gives for wavefunctions 
the relative deviation in real x-space less than < 10 -3 . 
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Needless to say that since the creation of quantum mechanics the one-dimensional quartic 
anharmonic oscillator 

d 2 

H = + rn 2 x 2 +gx 4 , (1) 

always attracted a lot of attention being among the most celebrated problem of quantum mechanics. 
The interest to these problems ranges from various branches of physics, from quantum field theory 
to chemistry and biology. It is especially true for the case when m 2 < and the potential has 
two minima. This problem is known in literature as the double-well potential. It was studied in 
hundreds papers, appeared practically in all books on quantum mechanics. A special emphasis was 
made to a domain m 2 — > oo - a domain where the barrier penetration is described by instantons 
(see e.g. 3,13]) that gives rise to the instanton physics. 

The first detailed study of (1) carried out by Bender- Wu at 1969-1973 revealed in this 
seemingly simple Hamiltonian the extremely rich analytic structure which looks intrinsic for any 
non-trivial eigenvalue problem of quantum mechanics and even though for quantum field theory. 
In fact, one of the most important unwritten conclusions was that in no way this problem can be 
solved exactly. The goal of the present talk is to give an approximate solution valid, actually, for 
any g > and m 2 . The solution is given in a form of the fairly simple expression for the ground 
state and the first excited state wavefunctions, which for any real x and for g > 0, real m 2 differs 
from the exact wavefunction for not more than a small number <5, 



^approximate ^ exact j 
^approximate 



< 5 . 



In our case the 5 ~ 10~ 6 . Evidently, it implies that any quantity related with the first two 
eigenstates like expectation values can be calculated with accuracy 5 2 . 

There are three basic analytic approaches to study the spectra in quantum mechanics: (i) per- 
turbation theory, (ii) WKB method and (iii) instanton calculus [lj]]. Each approach has its domain 
of applicability and usually these domains do not overlap. We attempt to combine (incorporate, 
unify) all three approaches into one by making interpolation. The most convenient object to in- 
corporate (i) and (ii) is the logarithmic derivative of the wavefunction. While the suitable object 
to incorporate the property (iii) is the wavefunction. A final form of the approximation depends 
on a few free parameters. Roughly speaking, their behavior as a function of g, m 2 is rather smooth 
and simple. They can be fixed variationally, although it is not very important: a small variation 
of the parameters does not lead to dramatic loss of accuracy. 
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As a first step to approach the problem let us remind the Symanzik rescaling for eigenvalues 
and eigenfunctions 



E(m 2 ,g) = g^E (^,2^j , m 2 , g) = * (xg 



1/6. 



g 2/3> 

It manifests that the original problem (1) is in fact one-parametric. The Hamiltonian ([I]) can be 
rewritten in the form 

d 2 

H = — r^r + ax 2 + 2x 4 , x G (-oo, +oo) , (2) 
car 

2 

where a = ^3 . This is the form of the Hamiltonian we are going to study with the real parameter 

a varying from negative to positive values. The Schroedinger equation for ([2]) reads 

d 2 ^f f +oc 
- — + ax 2 ^ + 2x% = E * , / \^\ 2 dx < 00 . (3) 

J — 00 

Eigenfunctions of ([3]) are sharply changing functions in x G R being characterized by a power-like 
behavior at |x| — ► and an exponentially-decaying one at \x\ — > 00. Following the oscillation 
(Sturm) theorem the nth eigenfunction has n simple (real) zeros. It seems natural to introduce 
the representation for eigenfunctions as follows 

= p n (x)e-^ , (4) 

where the phase <p(x) is a slow-changing smooth function and p n (x) is a polynomial of nth degree 
with real coefficients which has n real roots [eJ. Recently, it was obtained a remarkable result 6]: 
Any eigenfunction ^f(x) of (2) for any real a is entire function and it has infinitely-many simple 
complex zeros all situated on imaginary axis symmetrically! It implies that the phase <p(x) has 
infinitely- many logarithmic branch points in complex x-plane and has no singularities at real x. 
After substitution of Q into fl3J) we get the following equation 

y'_ y 2_ Pn-^P'n = E _ ax 2 _ ^ = , = (l og Vj/^))' , y(0) = . (5) 

Pn 

In order to define the problem ([5]) we impose two conditions that y has no simple poles at real x 
(i) and it grows at \x\ — > 00 not faster than polynomial (ii). The condition (i) implies that the 
coefficients of the polynomial p n (x) are those that the residues in the simple poles in the third term 
in l.h.s. of ([5]) vanish. The condition (ii) assures squire- integrability of the wavefunction It is 
evident that the polynomial p n (x) has parity p = (— V\ hence, can be written asp n (x) = x p P[«](x 2 ), 



and y(x) is odd, y{x) = y(—x) (for discussion, see 



a). 
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Prom the analysis of © it is easy to find asymptotic behavior of the phase 16], 

21/2 2, « , , / , 8E + a 2 1 fa A A 1 

if = —x \x\ + ^\x\ + (n + 1) log \x\ + 29/2 ^ + ^- - + . . . at |x| oo , (6) 

where ^4 is sum of squared of nodes, while 

E 2 £ 2 - a 4 2£(£ 2 - a) - 6 6 
^ = ~2 X + ~Y2~ X + 90 x + . . . at |x| - , (7) 

(see [§j]). It is important to note that the first two terms in ([6]) are defined by the equation ([5]) 
with omitted y' term, which is actually the Hamilton-Jacobi equation. Therefore, these two terms 
coincide with first two terms of the asymptotics of the classical action at |x| — ► oo. The third 
term in Q is also reproduced in the expansion of the classical action but with a wrong coefficient. 
The correct coefficient can be obtained if the first correction to the classical action is taken into 
account (quadratic fluctuations). The first three terms in JH]) grow when |x| tends to oo. In fact, 
they characterize a singularity at |x| = oo. These terms do not depend on energy E and are found 
explicitly. For the single well potential a > the expansion ([7]) is nothing but the (divergent) 
perturbation theory series near the minimum of the potential. It is an expansion around the true 
vacuum. 

Now let us construct a simplest function for phase (HJ) which interpolates small and large distance 

expansions, reproducing exactly the first three (growing) terms in ([6]). It has the form, 

A + (D 2 + 3a)x 2 + 4x 4 n + 1, n , , 

^ = 6(^ + 2sV 2 ^o g (D 2 + 2x 2 ). (8) 

The corresponding wavefunction dH) for n = (2k + p) excited state (for k = 0, 1, . . . and parity 
p = 0, 1) is equal to 

(M _ X PP k (x 2 ) j A + (D 2 + 3a)x 2 + 4x 4 | 

% ~ (D 2 + 2x 2 r^ eXP \ 6(D* + 2s»)i/» /' (J) 
(cf. [8] at k = 0), where A,D are free parameters and is a polynomial with real coefficients 
having positive roots only. There are two ways to find the polynomial either imposing the 
orthogonality conditions to the functions with smaller quantum numbers k = 0,1, . . . (k — 1), or 
requiring the absence of simple poles in the third term in l.h.s. of (|5|). Surprisingly, in concrete 
calculations these two conditions lead to polynomials whose coefficients coincide with high accuracy. 
Each function ^ is characterized by two free parameters A, D, which can be fixed if the function 
([9]) is taken as variational trial function. If k = in (|9|) both ground (p = 0) and the first excited 
(p = 1) states occur, respectively. Concrete calculations of the variational energies of these states 
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with ([9]) taken as a trial function lead to unprecedented accuracy [8[ (see below). Furthermore, if 
in the exponential in ([9]) in the term 4x (this term governs the asymptotic behavior of the wave 
function at large distances where the wave function is exponentially small) the factor 4 is replaced 
by a parameter, Wx 4 , minimization of such a trial function leads to a value of W which is equal 
to 4 with accuracy 10 ! 

In the case o < the vacuum at x = becomes the false vacuum and two classically degenerate 
vacua at x± = ±a 1//2 /2 appear. The expansion (J7J) becomes the expansion around false minimum 
(maximum) , it is not relevant physically. The expansion around one or another true vacuum should 
be considered instead, as a relevant one. This expansion can be easily derived and we skip it. 

In this case a new physical phenomenon of the quantum mechanical tunneling (barrier penetra- 
tion) occurs. There is a probability to meet the particle under the barrier, near x ~ 0. It decays 
exponentially when a — * —oo. This phenomenon is absent in ([9]). A prescription how to describe 
tunneling it is given in the celebrated Landau-Lifschitz book 9J]. It employs a linear superpositions 
of the wavefunctions centered at different minima which has positive parity for the ground state 
and similar one of negative parity for the first excited state. From the viewpoint of construction of 
the interpolation it is similar to making an interpolation between the expansion at one minimum 
and another one, and at |x| — > oo. Finally, such an interpolating function is a linear superposition 
of two off-centered functions Q, which can be written as 

_ Qttf) cogh « + x f A + + (D\ + 3a)x 2 + 4x 4 ] 

for the states of positive parity and 

_ C(^ 2 ) . <*-a f ^4- + + 3a)x 2 + 4x 4 ) 

% ~ (D 2 _ + 2x2)^+1 Smn (D 2 _ + 2x2)V2 eXp \ 6p 2 + 2x2)V2 J ' W 

for the states of negative parity. Parameter q + (q_) 'measures' a displacement of the peaks of 

wave function from the origin, x = 0. If a+ = the function (|10p becomes (|9"j). Here (q, ) is a 

polynomial of degree k with real coefficients having positive roots only. At fixed k any function 

depends on three free parameters a, A, D. 

In order to proceed further I need to remind two important and poorly known results: (i) a 

special form of perturbation theory in QM sometimes called 'Logarithmic Perturbation Theory' or 

'Non- linearization Method', and (ii) a connection between variational calculation and perturbation 

theory. 
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(i) A special form of perturbation theory is a certain iterative procedure developed for solving 
the Riccati equation ([5]) instead of the Schroedinger equation. For finding the wave function it is a 
multiplicative perturbation theory unlike a standard additive Rayleigh-Schroedinger perturbation 
theory. Such a multiplicative perturbation theory was developed for the first time by Price [l(| 
and then it was numerously rediscovered (for early history and discussion see [3] and references 

n 

therein). In presentation we follow closely to [TJ|. For simplicity we will consider the eigenstates 
for which the nodes are absent (ground state) or nodal positions are known. 

As a first step to develop the perturbation theory we make a choice of some square-integrable 
function an d calculate its logarithmic derivative 

y = (log* )' = ^ • (12) 
It is clear that is the exact eigenfunction of the Schroedinger operator with a potential 

Vo = ^ = y 2 -y'o, (13) 

where without a loss of generality we put their eigenvalue equals to zero, Eq = 0. It is nothing but 
a choice of the reference point for eigenvalues. Now we can construct a perturbation theory for 
Riccati equation taking an d yo, Vq as zero approximation, which characterizes the unperturbed 
problem. One can write the original potential V = m 2 x 2 + gx A as a sum, 

V = V + (V - Vq) = V Q + V X , (14) 

thus, taking a deviation of the original potential from the potential of the zero approximation as a 
perturbation. We always can insert a formal parameter A in front of V\ and develop a perturbation 
theory in powers of A, 

oo oo 

E = Y J ^ k E k , y = ^A fc y fc , (15) 

k=0 fc=0 

putting A = 1 afterwards. Perhaps, it is worth emphasizing that in spite of the fact that we study 
iteratively the equation ([5]), in general, this perturbation series has nothing to do with a standard 
WKB expansion. By substituting (|15|) into (0) we arrive at the equations which defines iteratively 
the corrections 

y'k ~ 2 yoVk = E k -Q k , (16) 

where 

Qi = V l , 
fc-i 

Qk = -^Vi- Vk-i > fc = 2, 3, . . . . 

i=l 
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It is interesting that the operator in the l.h.s. of (|16p does not depend on k, while Q k in the r.h.s. 
can be interpreted as a perturbation on the level k. The solution of ([To]) can be found explicitly 
and is given by 

Ek ~ r°° ,t,2 j — ' ( 17 ) 

J-oo ^ o dx 

Vk = % 2 r (E k -Q k )y 2 dx> . (18) 

J — oo 

It is easy to demonstrate that if the first correction y\ is bounded, 

\y\\< Const , (19) 

it provides a sufficient condition for this perturbation theory (|15p to be convergent Q]. Note that 
this condition is very rough and very likely can be strengthened. 

(ii) The first two terms in the expansion of energy (|15p in the above-described perturbation 
theory admit an interpretation in the framework of the variational calculus [5|] . Let us assume that 
our variational trial function ^q{x) is normalized to 1. We can calculate the potential Vq where 
^q{x) is the ground state eigenfunction and even put Eq = (see a discussion above). Formally, 
we construct the Hamiltonian Hq = p 2 + Vq for which Hq^q(x) = 0. The variational energy is 
equal to 

E var = j V'o-H'V'o = y i>oH ip + J ipo (H — Hq) ipp 

=Eo s v ' 

=Ei 

= E + E 1 (V 1 = V - Vq) > E exact . (20) 

Of course, ^q{x) could depend on free parameters. In this case both Vq and V\ depend on param- 
eters as well. Minimization of E var with respect to the parameters can be performed and the vari- 
ational principle guarantees that E var gives upper bound to the ground state energy. This simple 
interpretation (|20p reveals a fundamental difference between perturbation theory and variational 
calculus. Variational estimates can be obtained independently on the fact that the perturbation 
theory associated with trial function ^>q(x) is convergent or divergent. However, it seems natural 
to remove this difference by requiring a convergence of the perturbation series. In this case by cal- 
culating the next terms E2, E3, . . . in (|15p one can estimate the accuracy of variational calculation 
from one side and improve it iteratively from another side. An immediate criteria how to choose 
^o(x) in order to get a convergent perturbation theory is to have the perturbation potential V% to 
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be subordinate with respect to the non- vanishing potential of zero approximation Vq, 

< 1 , for |a?| > R . (21) 



11 

V 



From this point of view any function ([9j). (ll0p or, ([111) taken as an entry leads to a convergent 
perturbation theory. An open question is how to estimate the radius of convergency. 

The requirement (|2ip has a non-trivial physical implication: in order to guarantee a conver- 
gence of perturbation theory a domain where the wavefunction is exponentially small (classically- 
prohibited domain) should be reproduced as precise as possible. The same time a description of a 
domain where the wavefunction is of the order 1 is not important. It contradicts to a straightfor- 
ward physics intuition and underlying idea of variational calculus which, in particular, requires a 
precise description of the domain where the wavefunction is of the order 1. Needless to say that 
namely the latter domain gives a dominant contribution to the integrals which define the energy 



121. 



in the variational calculations. Similar conclusion was presented in 

Results. 

Ground state 



a = l 

This is the case of a single- well potential (anharmonic oscillator). The variational parameters in 
(fTUl) are 

D = 4.33441 , A = -9.23456 , a = 2.74573 . 

Ground state energy is 

Evar (= ^o + ^i) = 1.607541302594, 
while the first correction to it is 

AE var (= E 2 ) = -1.2552 x 10~ 10 . 

Eventually, 

E V ar = E var + AE var = 1.607541302469 , (22) 

where all 13 digits are correct, since the next correction E3 ~ 10 -14 . The rate of convergency 
seems extremely high, ~ 10~ 4 ! The dependence the energy on a is very weak: the results are 
almost unchanged if a = (and ([TO]) becomes ([9])). It is the case for a > 0. 
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FIG. 1: Logarithmic derivative yo = (f' int (see Q) as function of x for double- well potential (2) 
with a = —1 



a = -1 

This is already the case of a double-well potential. The variational parameters in (|10p are 

D = 4.059888 , A = -12.4816 , a = 3.07041 . 

In comparison with the case a = 1 the parameters D, a are slightly changed unlike the parameter 
A which decreased in ~ 40%. Ground state energy is 

E var (= Eq + Ei) = 1.029560832093, 

while the first correction to it is 

AE var (= E 2 ) = -1.0382 x 10~ 9 . 

Eventually, 

E V ar = E var + AE var = 1.029560831054 , (23) 

where all 13 digits are correct, since the next correction E% ~ 10~ 13 . Similar to the case a = 1 the 
rate of convergency seems extremely high, ~ 10 ! 

On Fig.l one can see the behavior of logarithmic derivative yo vs x > 0. It is a very smooth 
function. The first correction \yi\ (see Fig. 2) has very interesting behavior: it is of the order of 
10 -4 at x < 1 (in the domain which give a dominant contribution to the energy integral (|20p ). then 
it starts to grow and reaches the maximum ~ 0.006 at x ~ 3.9 in the domain which gives negligibly 
small contribution to the energy integral, yo(x = 3.9) ~ 15.2 and "ifo(x = 3.9) ~ 10~ 9 . In this case 
5 ~ max\y\\ ~ 0.006 (see p.l). Notably, y\ oc 1/x 2 at x S> 1. Similar behavior is demonstrated 
by the higher corrections = 2,3: they are very small at x < 1, its maximum is reached for 

x > 1 but its position is systematically reduced with n. 
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FIG. 2: The first correction y\ for a = —1 (cf. Fig.l) 

(X — &crit 

Let us consider the Hamiltonian (1). The ground state function is symmetric w.r.t. x — > —x. 
Hence, it has to have an extremum at x = 0. For any fixed g there exists a value m 2 rit < such 
that for to 2 > m 2 rit this extremum is a maximum, otherwise a minimum. It is easy to find out 
that the critical point m 2 rit corresponds to the vanishing ground state energy, E = [l?]]. Using 
the function (fT0|) it was calculated the critical value m 2 rit = —2.2195970861 for g = 1. For the 
Hamiltonian ([2]) 

E(a crit = -3.523390749) = . 

It is quite interesting from physical point of view that for a family of double-well potentials with 
fixed g there exists a domain > m 2 > (m 2 ) cr it where the ground-state eigenfunction has the 
maximum at the origin, which corresponds to the position of the unstable equilibrium similar to 
what takes place for the single-well case. It implies that the particle in such a potential with the 
ground state energy above the barrier, E > 0, somehow does not feel the existence of two minima. 
In this domain WKB consideration at x ~ is not valid. 

a = -20 

This is the case of a double-well potential. The variational parameters in (fT0|) are 

D = 6.765663 , A = -286.6456 , a = 49.6136 . 

Ground state energy is 

En* (= £ + £i) = -43.7793127 , 
while the first correction to it is 

AE var (= E 2 ) = -3.81 x 10~ 6 . 
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Eventually, 

E var = E var + AE var = -43.7793165 , (24) 

where all 9 digits are correct, since the next correction E% ~ 10 -8 . Similar to the cases a = ±1 the 
rate of convergency seems still high, ~ 10 -2 ! 

First excited state 

a = -20 

This is the case of a double-well potential. The variational parameters in (jllh are 

D = 5.584376 , A = -246.64375 , a = 38.82768 . 
The parameters seems quite close to those for a = 20 for the ground state. Ground state energy is 

E var (= Eq + Ex) = -43.77931637 , 
while the first correction to it is 

AE var (= E 2 ) = -9.3618 x 10~ 8 . 

Eventually, 

E var = E var + AEvar = -43.77931646 , (25) 

where all 10 digits are correct, since the next correction E% ~ 10~ 10 . Similar to the cases a = 20 
for the ground state the rate of convergency seems still high, ~ 10 -2 ! 

Energy Gap 

By definition the energy gap is 

AE — Eji rs i excited state Eg rouna state • 



otential (2) at a — ► — oo it can be calculated and in one-instanton approxi- 



For the double-wel 
mation it reads 

2 11/4 , ,c, A V2M 3/2 / 71 1 6299 1 2691107 1 2125346615 1 



AE 



— F a 5, *e s 1- 



y/if ] " ] " \" 12 V2|a| 3 / 2 288 2|a| 3 10368 2 v / 2|a| 9 / 2 497664 4|a| 6 " ' 

(26) 



The first term in the expansion is a common knowledge. All other coefficients in the expansion 
are due to J Zinn- Justin, 1981-2005 (see [3| and references therein), they were obtained using the 
so called exact Bohr-Sommerfeld quantization condition. The coefficient 71/12 was independently 
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14l ] and references 



calculated (and confirmed) in two-loop instanton calculation by E Shuryak (see 
therein). It seems highly desirable to perform three-loop calculation to check the next coefficient 
in the expansion. It is an asymptotic expansion 13]. 

It is interesting to compare the energy gap AE at a = —20 calculated in the (convergent) 
perturbation theory with (|10p and (jlip as zero approximations and with use of the (asymptotic) 
expansion (126j) . Subsequent expressions show how AE evolves from pure variational results to ones 
with the first corrections E2 taken into account and then to ones with the second corrections £3 
involved 

AE var = 1.03282 x 10~ 7 , 
AE^ r = 1.06529 x 10~ 7 , 

AE$ r = 1.06525 x 10~ 7 . (27) 
In the last expression AEi,a r all six significant digits are correct. Now how the energy gap AE 



looks like as function of a number of corrections to one- instanton result included (see (|26p ) 

one — instanton = 1.12154 x 10 _ (5.3% deviation) 

one — instanton + 1st correction = 1.06908 x 10~ 7 (0.36% deviation) 

one — instanton + 1st and 2nd corrections = 1.06754 x 10~ 7 (0.22% deviation) 

one — instanton + four corrections = 1.06738 x 10~ 7 (0.20% deviation) (28) 

The numbers in brackets are relative deviations from (|27p . The fact that a deviation stays almost 
the same after adding 2nd, 3rd and 4th corrections likely indicates the maximal accuracy based on 
a use of asymptotic expansion is reached and in any moment the result can blow up. A comparison 
of two values (|27p and (|28p shows that they do not agree in the 4th digit. 

In a conclusion I have to say that a similar consideration based on interpolation of phase between 
small and large distances was done for sextic oscillator and the Zeeman effect on hydrogen. In both 
cases the exceptionally high accuracies were obtained. 
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